Waveform- Controlled Terahertz Radiation from the Air Filament Produced by 

Few-Cycle Laser Pulses 
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Waveform- controhed Terahertz (THz) radiation is of great importance due to its potential apph- 
cation in THz sensing and coherent control of quantum systems. We demonstrated a novel scheme 
to generate waveform- controlled THz radiation from air plasma produced when carrier-envelope- 
phase (CEP) stabilized few-cycle laser pulses undergo filamentation in ambient air. We launched 
CEP-stabilized 10 fs-long (^ 1.7 optical cycles) laser pulses at 1.8 /xm into air and found that the 
generated THz waveform can be controlled by varying the filament length and the CEP of driv- 
ing laser pulses. Calculations using the photocurrent model and including the propagation effects 
well reproduce the experimental results, and the origins of various phase shifts in the filament are 
elucidated. 

PACS numbers: 52.38.Hb, 41.60.-m, 42.65.Re, 32.80.Fb 



Air-plasma based terahertz (THz) wave generation at- 
tracts much attention as it provides unique tools for non- 
linear spectroscopy, imaging and remote sensing [l|-(3|. 
As an intense femtosecond laser pulse undergoes filamen- 
tation in ambient air, axial polarized THz radiation in 
forward direction was generated due to the transition- 
Cherenkov-type radiation [4]. More intense THz radia- 
tion field (up to MV/cm range) can be produced by us- 
ing two-color laser field, where four wave mixing model 
and photocurrent model are used to explain the underly- 
ing physics . The spatiotemporal dynamics of THz 
emission from the air plasma generated by the two-color 
laser field has recently been discussed and the scheme 
of tailoring the THz emission spectrum was proposed by 
adjusting the tunneling ionization events P. [iq|. 

Alternatively, intense THz radiation can be generated 
from air plasma driven by an intense few-cycle laser pulse. 
By measuring the amplitude and polarity of generated 
THz emission, the carrier-envelope-phase (CEP) of few- 
cycle laser pulses can be inferred [11]. The THz emission 
was interpreted as a result of asymmetric ionization in- 
duced quasi-dc current produced by the few-cycle laser 
pulse, as confirmed in subsequent theoretical investiga- 
tions [13, [13. 

It is well understood that a converging Gaussian beam 
experiences a phase shift of tt, the so-called Gouy phase 
shift through the laser focus, which has been measured 
using the stereo above-threshold ionization (ATI) method 
in high vacuum conditions [M]. However, the existence 
of Gouy phase shift in femtosecond laser filaments in air 



remains controversial 15 



Recent investigations show that waveform-controlled 
THz radiation is of great importance due to its poten- 
tial application in THz imaging and coherent control of 
molecular dynamics 17, 18|. Can we control the wave- 



form of femtosecond laser filament based THz radiation? 
Understanding the laser phase effect and propagation ef- 
fects on the THz radiation from the few-cycle laser pulse 
produced filament is crucial for the THz generation based 
CEP metrology jTl| and the generation of intense and 
waveform-controlled THz radiation. 

In this Letter we propose and demonstrate a novel 
scheme to generate waveform- controlled THz radiation 
from the air filament produced by few-cycle laser pulses. 
Variation of THz waveform and even its polarity inver- 
sion are found in the spatially-resolved measurement of 
THz emission along the filament. THz waveform can be 
controlled by varying the filament length and the initial 
CEP of the driving laser pulses. 

The CEP stabilized infrared (IR) few-cycle laser pulses 
are produced using a home-built optical parametric am- 
plifier (OPA) system and a hollow fiber compressor, 
schematically shown in Fig.lTJa). The OPA and compres- 
sor are in principle similar to our previous system [l9| 
except that the pumping laser is the Tiisapphire laser 
amplifier facility (Coherent Elite-HP-USX) with pulse 
energy of 5.2 mJ, duration of 25 fs and repetition rate 
of 1 kHz at 800 nm. The major portion of the output 
beam, 4.6 mJ, is used to pump the OPA system while 
the left as a probe pulse for sampling the generated THz 
temporal waveforms. In the three-stage OPA amplifier, 
the idler beam of the 3^^ stage OPA amplifier is used 
as the output, so its CEP is passively stabilized with a 
rms jitter of ~ 400 mrad. The 40 fs-long CEP-stabilized 
laser pulses at 1.8 fim are then compressed down to 10 
fs (~ 1.7 cycles) with the maximum pulse energy of 450 
/iJ. When the laser pulse is focused into ambient air by 
using a spherical mirror f = 150 mm, a stable lumines- 
cence filament of ~ 12 mm- long is formed by using 300 
/iJ pulse energy selected by an iris diaphragm, as shown 
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FIG. 1: (a) Schematics of the experimental layout; (b) Image 
of the filament formed by focusing the few-cycle laser pulse 
into ambient air; (c) Schematics of the filament length control. 



in Fig. [TJb). The generated THz radiation are measured 
by using the standard balanced diode geometry electric- 
optic (EO) sampling technique. 

In order to investigate the evolution of driving laser 
pulse and THz radiation in the filament, we insert a sharp 
stainless steel blade (~ 0.2 mm x 4 mm x 20 mm) into 
the plasma column to stop the filament, which can be 
moved by a motor stage along the laser propagation di- 
rection defined as the z coordinate henceforth. As shown 
in Fig. m^c), the blade is positioned to allow the laser 
beam hit onto its sharp edge, so that the detection of 
THz radiation from the left filament is lest influenced. 
In so doing, the integrated THz signals from filaments of 
different lengths are measured. We calibrate experimen- 
tally the deviation in the collection efficiency of the THz 
signal detection optics due to the change in the distance 
between the filament and the detector. 

Fig.[2fa) shows two THz waveforms by blocking the fil- 
ament at the position of 3 mm (the distance from the vis- 
ible starting position of the filament to the blade) and 10 
mm, respectively. It is reasonable that the measured THz 
amplitude from the 3 mm-long filament is much smaller 
than from the 10 mm-long one, since the measured signal 
comes from the integrated emission from the unblocked 
filament. One can also note that the recorded THz radi- 
ation from the 10 mm-long filament reverses its polarity 
comparing to that from the 3 mm-long filament, indicat- 
ing the phase variation of the driving laser field in the 
filament. By moving the blocking blade continuously in 
a step length of 0.5 mm, the THz waveforms as a function 
of filament length are recorded and plotted in Fig. [2fb). 
One can see that, the THz emission gradually changes 
its amplitude and polarity when the length of filament is 
varied, and at certain length (~ 5 mm-long filament) the 
THz emission signal goes to zero and then reverses the 
polarity thereafter. One can therefore control the THz 
waveform by varying the filament length. 

For analyzing the amplitude modulation of THz radi- 
ation, the amplitude values when the THz signal is the 
maximum, at ~ 1 ps in Fig.[2jb), are plotted in Fig.[2jc). 
With the help of the above-mentioned calibration of col- 



lection efficiency, the spatially resolved THz emission as 
a function of the position in the filament can be retrieved, 
also shown in Fig. EJc). One can see that the THz emis- 
sion first increases in one (negative) polarity and then 
decreases to zero at ^ 3 mm, after which the amplitude 
increases in the other (positive) polarity and decreases at 
the tail of the filament. 

The polarization of THz emission was measured to be 
nearly linear so that we excludes that the generated THz 
emission is originated from transition- Cherenkov- type ra- 
diation. We consider the THz emission originated from 
the transient photocurrent driven by the propagating in- 
tense few-cycle laser fields in the plasma [6|. The prop- 
agation of intense few-cycle laser fields in a dispersive 
medium can be described by using a propagation equa- 
tion in an axial-symmetric coordinates [20|, 



d^E{r, z, uj) = [—- vi ^ik(yj)\E{r, z, uj) 
-Pnl — — T^^Ji( 



(1) 



where E{r^z^uj) is the frequency domain laser field, the 
nonlinear polarization Pnl accounts for the Kerr effect 
and the polarization Jioni is caused by photoelectrons 
from the tunneling ionization of N2 and O2 in air. The 
collective motion of the tunneling ionized electrons re- 
sults in a directional nonlinear photocurrent surge de- 
scribed by 



dtJe{r,z,t) + VeJe{r,z,t) = —pe{r,z,t)E{r,z,t) (2) 

m 
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FIG. 2: (a) THz temporal waveforms measured from the fil- 
aments of different length 3 mm (red solid line) and 10 mm 
(black dashed line), respectively; (b) THz waveforms mea- 
sured by moving the block continuously along the filament; 
(c) The measured THz amplitudes (green solid squares) taken 
from Fig. 2(b) at the delay time of 1 ps, the calculated results 
(black solid line), and the retrieved THz amplitude (red solid 
circles) as a function of the position in the filament. 
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FIG. 3: (a) The calculated THz amplitude (red solid line) 
and asymmetry parameter (blue dashed line), X)(J^ — J^), 
in air plasma filament; (b-d) At positions of A, B and C, 
the laser field envelopes (blue thin line), carrier oscillations 
(black thick line), electron densities (green dashed line) and 
transient photocurrents (red dotted solid line). 



where Ve^e^m and pe denote the electron- ion collision 
rate, electron charge, mass, and electron density, respec- 
tively. The transient current at each propagation step of 
the calculation is treated as radiation source of far field 
THz emission (2lj 



An Co J c^R 



(3) 



where R is the distance between the THz point source 
and the detection plane. Eq. (j3j) describes the THz ra- 
diation towards all directions. The forward far field THz 
emission near the propagation axis from the unblocked 
air plasma has been added up to simulate the measured 
results. The initial few-cycle laser field is of a Gaus- 
sian beam with the waist size of wq = 260 /im (1/e^) at 
center wavelength of 1.8 /im, duration of r = 10 fs in 
full width at half maximum (FWHM) and pulse energy 
of 300 /iJ. By optimizing the initial CEP of the driving 
few-cycle laser field to be (/:5o = 0.33 tt, the THz amplitude 
as a function of filament length is found in consistence 
with the measured results, as shown in black solid line in 
Fig. He). 

For a further understanding of the THz emission char- 
acteristics in the air plasma, we show in Fig. [3ja) the 
space-resolved THz emission calculated from Eq. (|3|) 
along the filament. As one can see that the THz ampli- 
tude modulates from negative polarity at the beginning 
of plasma to positive one at the end, which results in the 
inversed THz emission polarity shown in Fig. [2fc). In a 
simplified picture, the d.c. transverse current is related 
to discrete electron ionization events and the pulse vector 



potential A{ti) modulation 13, [22^, which is described by 
Je cx 6piA{ti). In order to distinguish electrons with 
opposite drifting velocities, we define two quantities J+ 
and J~, which describe the positive current density and 
negative current density, respectively. We calculate the 
asymmetry of the summed current densities "^{J^ — Jq) 
which matches the THz emission as shown in Fig. [3fa). 

The laser fields, electron densities and transient pho- 
tocurrents at three positions labeled as A, B and C in 
Fig. [3l^a) are calculated and plotted in Figs. [3l^b)- [3l^d), 
respectively. At the position of A, where the THz emis- 
sion is the maximum in the negative polarity, the ampli- 
tude of negatively current density J~ is larger than the 
positively current density J+ following the laser field os- 
cillation. At the position of B, the negative current den- 
sity is less than the positive one so that the THz radiation 
is positively polarized. At the position of C, the almost 
equal current densities in the oscillation result in the near 
zero THz radiation. From Figs. [3l^b)- [31(d), one can also 
note that the increments of produced electron density pe 
are mainly from two attosecond bursts in the opposite 
directions, that are originated from the tunneling ioniza- 
tion in a single optical cycle, making the THz emission 
extremely sensitive to the variation of driving laser fields. 
The attosecond electron bursts can be easily controlled 
by varying the intense few-cycle fields, comparing to the 
THz radiation generation by two-color multi-cycle laser 
fields [9]. 

Since the THz emission is shown to be sensitively de- 
termined by the intense few-cycle laser field, we then in- 
vestigate the CEP variation in the air plasma. As in- 
dicated previously [23], the CEP is the difference of the 
pulse phase, i.e. carrier phase, and the pulse front. The 
carrier phase shift in a dispersive medium can be written 
as 

(l^uoir^z) = n{r,z)^z^(j)G{z,Zeff)^(l)off-a{r,z) (4) 

where the first term is the phase shift at the carrier 
frequency induced by spatial dependent refractive index 
n(r, z), the second term is the Gouy phase shift with the 
effective focal length ^e//, and the third term denotes 
the off- axial phase shift due to diffraction effects. For 
simplicity we only consider the axial phase shift that in- 
cludes the first two terms. We define the second term 
as the modified Gouy phase shift that is characterized 
by Zeff^ since the Rayleigh length is no longer valid in 
filament. 

The calculated modified Gouy phase shift decreases 
with a relatively slower slope than the Gouy phase shift 
under the linear focus condition, as shown in Fig. lU^a). 
The modified Gouy phase shift due to the extension of 
filament confirms the previous theoretical prediction that 
the Gouy phase shift stems from transverse spatial con- 
finement of a finite beam [2^. The CEP of laser fields 



is also determined by the propagating pulse front, which 
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FIG. 4: (a) The modified Gouy phase shift (red dashed line) 
and the Gouy phase shift (black solid line) of linear focus con- 
dition; (b) The pulse front shift (blue solid line) and the CEP 
shift (red dotted solid line) of the laser field in air plasma. 



is defined as Eq. (17) of the Ref. [23^ and can be re- 
trieved from the calculation results. The axial pulse front 
is shown in Fig. IH^b). As one can see the pulse front 
experiences a jump around the center of the plasma re- 
gion (z - -7 mm ~ 4 mm), as a result of the plasma- 
induced dispersion. The total axial CEP shift is ~ 1.4 
TT and possessing a hump structure in the plasma region, 
which is originated from the modified Gouy phase shift 
and plasma effects. 

It should be noted that the modified Gouy phase shift 
also applies in the filament produced by multi-cycle laser 
pulses, in which the CEP is not a significant parameter 
for THz generation. In the recent studies on THz Air- 
Biased-Coherent-Detection (ABCD) |2f 



26|, the Gouy 



phase shift of the THz emission and the weak ultrafast 
probe pulses around the focus was observed. The present 
work reveals that, due to the nonlinear interaction in air 
plasma, the intense laser field experiences not only the 
modified Gouy phase shift but also the variation of pulse 
front. 

Furthermore, we look into the initial CEP, (po^ depen- 
dence of the THz emission from the air plasma, cpo is 
adjusted by passing the laser beam through a pair of thin 
wedges mounted on a motor stage. By increasing the ini- 
tial CEP gradually in a step size of 0.27r, we obtained 
a series of 2-dimentional THz waveform maps similar to 
that shown in Fig.[2fb). The amplitudes of THz emission 
along the filament are retrieved and plotted as a function 
of (fo in Fig. [5ja). We found that the THz waveforms 
change back to be the same as those shown in Fig. [2lb) 
after a change of initial CEP by 27r. This indicates that 
in air plasma the propagation effect on the CEP vari- 
ation is independent of the initial CEP of driving laser 
fields, which is also verified by our numerical simulations. 
The phases of few-cycle pulses experience a fixed shift of 
over TT, which is consistent with the CEP shift through 
the air plasma. This observation helps to validate the 
CEP metrology based on THz generation [ll|. How- 



ever, due to the variation of CEP through the plasma, 
an offset of CEP must be determined and taken into ac- 
count for measuring the CEP of laser pulses. As shown 
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FIG. 5: (a) Measured and (b) simulated THz amplitude as 
functions of relative CEP shift and filament length. 



in Fig.[5fb), the simulated results using the photocurrent 
model and propagation equation agree well with experi- 
mental results. Since the polarity of THz radiation can 
be controlled by varying the initial CEP of driving laser 
pulses, one can maximize the output THz radiation in- 
tensity from the filament by optimizing the initial CEP 
to minimize the intensity cancellation due to opposite 
polarities. 

In conclusion, we have demonstrated a scheme to gen- 
erate waveform-controlled THz radiation from air plasma 
produced by CEP-stabilized few-cycle laser pulses. The 
generated THz waveform can be controlled by varying the 
length of the filament and the CEP of the driving laser 
pulses. The waveform evolution of THz radiation in the 
filament is due to the characteristic carrier phase shift 
and variation of laser envelope during propagation. Due 
to the ability to control the spatial location of femtosec- 
ond laser filamentation using pre-chirped laser pulses, the 
demonstrated scheme of waveform-controlled THz radi- 
ation generation is suitable for remote sensing exploring 
the phase property of THz radiation. The knowledge 
about the phase shift in a filament would also be valuable 
for optimizing the generation of phase-matched attosec- 
ond pulses driven by intense few-cycle laser pulses. 
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